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Abstract 



We consider the branching ratio and the CP asymmetries in B — > fQ(980)Ks decay to the end of deter- 
mining the deviation of the time-dependent CP asymmetry from sin(2/3), ASf^^Xs = "Vo^s'^/o^s "^^(2/^)' 
arising from Standard Model physics. We obtain ASj^Ks within the context of the QCD factorization 
framework for the B — > fQ{980)Ks decay amphtudes assuming the /o(980) is a qq state and employing a 
random scan over the theoretical parameter space to assess the possible range in ASf^^g. Imposing the 
value of the experimental branching ratio within la and 3a, respectively, of its central value as a constraint, 
we find the range of ASf^Xs [0.018, 0.033] for a scan in which the parameters are allowed to vary 

within la of their central values and the range [—0.019, 0.064] for a scan in which the parameters vary 
within 3a of their central values. 
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I. INTRODUCTION 



In the Standard Model (SM), all CP- violating effects derive from a single, complex phase of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix and predicate a distinctive pattern of CP-violation [Tj . 
For example, in the decay of a B-meson to a CP-eigenstate /, the time-dependent asymmetry Sf 
realized from b — > sec decay determines sin(2/5), where is given by = arg(— V^^V2,/V^d^b) and 
Vij is a CKM matrix element [21 El H]. This quark-level transition can be studied in a variety of 
B-meson decays, and departures of the determined time-dependent asymmetry from sin(2/5) could 
signal the presence of non-SM physics, which may occur m. B — B mixing, in the decay amplitude, 
or in both [5]. 

In this paper, we consider the decay B —>■ fo{980)Ks, which is mediated by the b —>■ sqq transition 
at one-loop-order in the weak interaction. The decay B — /o(980)-ft'5' is one of several penguin- 
dominated modes which probe sin(2/?). In contrast, in i? — >■ J/ipKs decay, and related charmonium 
modes, the b — > sec transition operates at tree level. Were the time-dependent asymmetries in 
tree- and penguin-dominated modes to differ, then non-SM physics could be at work in the penguin 
process [H]. Current experimental results suggest that this could be the case [7j, though definite 
conclusions require both experimental results of improved precision and theoretical estimates of the 
subleading SM corrections. The numerical size of the SM corrections depend on the specific decay 
mode, mimicking the appearance of non-SM physics [8], so that the needed estimates demand some 
care. In this context it is worth noting that the "wrong phase" penguin contribution, proportional 
to the weak phase of b ^ suu, is particularly small in B ^ J/ipKs decay; indeed, the deviation of 
the time-dependent asymmetry Sj/^Ks from sin(2/3) is C(10^^) [HI |10] — it is suppressed by both 
CKM and loop effects. Thus the comparison of this asymmetry to a "tree-only" determination of 
sin(2/5) permits a sensitive assay of new physics in B^ — B^ mixing. Currently this last is consistent 
with Sj/^Ksi as well as with other determinations of sin(2/5) which employ information on the 
sides of the unitarity triangle, at the 0(10%) level [HI [12]. In the case of the penguin modes, the 
wrong-phase penguin is larger as it is suffers only CKM, i.e., 0{\^) ~ 0.04, suppression. In these 
modes the computed SM deviations from sin(2/3) determine a much-needed baseline against which 
the experimental results can be assessed for new-physics effects, as new-physics-induced deviations 
from sin2/3 could certainly be channel-dependent as well. Systematic studies of the SM corrections 
exist [ini HH US] in a variety of modes. The B /o(980)i^5 mode has received less attention, 
perhaps due to the ill-known quark structure of the /o(980) [161 E]j this is a deficiency we wish to 
remedy. 

The /o(980) is a fairly narrow resonance of non-Breit-Wigner form which couples to vrvr and KK 
final states. The quark structure of the /o(980) meson is not well established. Much discussion has 
revolved around whether it is better regarded as a q^(f state [TB] or, perhaps, as a KK molecule [IHl 
[20] . The theoretical correction relevant to the interpretation of Sj^Ks as a measurement of sin2/3 
comes from the presence of the b —>■ suu transition in the decay amplitude, and interconnected 
issues complicate its analysis. Not only must we consider the possible non-gg structure of the /o 
resonance, but we must also recall that strong final-state interactions exist in vrvr scattering in 
the J = 0, J = channel [211 ESI ES]. Such final-state interactions give rise to the scalar form 
factor, and particularly the partial width r(/o(980) nir ,KK), which can be computed through 
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the unitarization of a scattering kernel compatible with low-energy constraints [2^ [25l 126] . We 
presume the /o(980) resonance to be sufficiently narrow relative to the energy released in the B 
decay that we can approximate the full decay amplitude as the product of the two-body decay 
amplitude A{B fo{980)Ks) with r(/o(980) ^ tttt^KK), where we refer to Ref. [21] for a 
discussion of the assumptions implicit to this treatment. In particular, the strong phases associated 
with the long-distance physics of /o decay are presumed to be universal and not modified by the 
B-decay environment, so that the decay-specific phases are captured by the application of QCD 
perturbation theory in the heavy quark limit to the two-body decay process. Corrections to the 
picture we employ can be estimated, though not in a systematically improvable way. For example, 
so-called final-state rescattering has been estimated for various b ^ s penguin modes and 
OZI- violating effects can also be considered, the latter contributing significantly to J ftp —>■ (pTrn 
decay, e.g. [251 126]. We reserve discussion of the impact of the computed /o nn and fo—*KK 
partial widths, as well as of possible OZI- violating effects, to a subsequent publication [2H|- Our 
particular focus is in the study of the sensitivity of our prediction of the CP asymmetries to the 
assumed quark structure of the /o(980); the latter enters in the evaluation of the hadronic matrix 
elements of the b — > suu transition. A priori one might think such effects to be most important 
in the assessment of the deviation of Sj^Ks from sin(2/?). We proceed in the same vein as Cheng 
and collaborators [161 IE] 5 treating the /o(980) as a qq state and employing the QCD factorization 
approach [291 EO] for the hadronic matrix elements. We differ from this earlier work in the treatment 
of the i? — *• /o form factor. We then perform a random scan over the theoretical parameter space, 
after Beneke [T3|, to assess, in part, the sensitivity of the results to the employed hadronic matrix 
elements and hence to the implicitly assumed structure of the /o(980). 

Let us give a brief outline of this paper. We begin, in Sec. [n| with a description of the theoretical 
framework and briefly review the role of the b — > sqq transition, with q E u,d, in the determination 
of the CP- violating parameters Sf and Cj, as well as of the two-body branching ratio. We present 
a synopsis of the pertinent QCD factorization formulae as well. In Sec. Ill we present the input pa- 
rameters we use in our numerical calculations, reporting the CP-violating parameters which result 
from the use of our defined "default" set of input parameters. It should be emphasized that the fac- 
torization theorem from which the QCD factorization approach follows holds only at leading power 
in the heavy quark mass, and the estimate of l/rrib power suppressed terms, though of apparent 
phenomenological importance, is uncertain. Thus our random scan over the theoretical parameter 
space, effected to explore the possible range of ASj^Ks ^foKs^ takes both uncertainties in the 
theoretical inputs and in the assessment of the 0{Aqcd/Mb) corrections into account. We report 
these results in Sec. Ill as well. We conclude with a summary of our results and an outlook on 
future work in Sec. [IVl 



II. THEORETICAL FRAMEWORK 



We pattern our analysis after that of Beneke, Ref. [H], and, indeed, adopt a common notation. 
In particular, we employ the QCD factorization framework [291 EO] for the computation of the 
B — > foKs matrix elements [17] and perform a random scan over the space of possible input 
parameters to estimate the uncertainty in their computation. We begin by recalling that the CP 
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asymmetry Af{t) into a CP eigenstate / is given by 



Here we have neglected AF where AF = Th — is the width difference of the B eigenstates. 
We note AMb = Mh — Ml is the mass difference of the B eigenstates, 5/ is the CP asymmetry 
generated by the interference oi B - B mixing and direct decay, and C/ is an asymmetry reflective 
of direct CP violation. We recall pQ 

where 

'q\ AjB^ f) 
^pJ^AiBO^f)^ 

with the factor {q/p)B characterizing B — B mixing. If we treat the /o(980) resonance as if it were 
a stable particle, as in Ref. [T7], then / = /o(980)ii's, and we write 



^foKs 



where the factor {q/p)K characterizes K ~ K mixing. In this event, / is a two-body final state, and 
we write the decay amplitude as 

A{B ^ /) = Ae a^^ + A„ oc (1 + e"^^ df) , (5) 

with df = |Au/Ac|(a.j/a.j) and Xq = VqbV*^ for g G m, c to determine 

^ ^ ^ ^' l + 2Re(rf/)cos7 + ' ^' 

where f]f is the eigenvalue of the CP-operator associated with the eigenstate /. Moreover, we have 

^ 2lm{df)sm-f 
^ l + 2Re(rf/)cos7+ |c//|2 ' ^ ' 

so that if \df \ is small, the functions ASf and Cf show little correlation [Hj. The value of ASf 
in the SM shows certain systematic trends with / For example, if a color-suppressed tree 

amplitude C contributes to a^, then ASf is much larger than it would be if it were absent. In the 
latter event, the w-quark penguin amplitude drives aj. Generally, df oc (±C -|- P'^)/P^. If the 
parameters describing the B ^ f decay amplitude make \P^\ small relative to \C\, then ASf can 
range over a wide array of values. Such large excursions can be controlled, however, by demanding 
that the amplitudes be consistent with the empirical branching ratios [14]. In the case of current 
interest, B foKs decay, there is no C amplitude, so that we expect ASfgXs to be small on general 
grounds, and the imposition of a branching ratio constraint should no longer be crucial. However, 
we have found exceptional regions in theoretical parameter space for which P'^ is small with respect 



4 



to and thus ASj^Ks large, so that it is, in fact, crucial to apply the branching ratio constraint 
to eliminate these large excursions. 



If we treat 



/o (980) -ft''' decay as a two-body process, then the decay rate is 

P 



8ttM^ 



where M = 



foK^), as we shall describe in detail, and 



P 



2M, 



(9) 



B 



We recall that the branching ratio is given hj T/Tb, where is the total decay width of B meson. 
In reality, the /o(980) is not a stable particle; rather, it is a resonance which decays to both vrvr and 
KK final states. We wish to investigate the role of finite-width effects explicitly in a subsequent 
publication ^28|: we neglect them here. To begin, we rewrite the two-body decay amplitude as 



A{B ^ foKo) = A^^X,, + A^X, + A^X^ + A^X, 



(10) 



where the superscript (n, s) refers to the non-strange and strange quark components of the /o(980) 
respectively, so that 

A 



d 



foKs 



Ac 



A- + AI 



11) 



We now proceed to calculate the A^'^ amplitudes. 



A. B ^ JqK Decay in QCD Factorization 



The decay amplitudes of exclusive hadronic S-meson decays can be systematically analyzed in a 
combined expansion of inverse powers of the heavy quark mass rrib and the strong coupling constant 
as [311 E21 |33]. The QCD factorization approach, in specific, permits the rigorous computation 
of these amplitudes, for certain two-body final states, in leading power in 0{AQCD/^b) and in a 
power series in as{fi), where /i ~ 0{mh) ^29]. Its starting point is the effective weak Hamiltonian 
for charmless hadronic B decay, consisting of a sum of the products of CKM matrix elements, local 
operators Qi, and Wilson coefficients Ci{^), evaluated in next-to-leading order (NLO) precision 
in asifi) [3l]. Some dissension in the literature exists concerning the ingredients of the leading 
power analysis [55] , particularly in regards to the charm-quark penguin contributions [5^1 EZ] : in 
the QCD factorization approach any non-factorizable charm-quark penguin contributions appear 
in 0{AQCD/^b) corrections. Although recent work suggests that such charm-quark penguin effects 
may be needed to explain the -B-decay data to light charmless mesons [3Hj, note also Ref. [SH], in 
our context we can safely neglect such developments, as they cannot make ASf^Xs larger. 5-meson 
decays to scalar- and pseudoscalar-meson final states have recently been studied by Cheng and 
collaborators in a series of papers [151 [HI [Hlj they employ the QCD factorization approach and 
treat the /o(980) resonance as a qq state. Their results connect to those of Beneke and Neubert for 
-B-meson decays to pseudoscalar- and vector- meson final states [30j, as the amplitudes for the scalar- 
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and pseudoscalar-meson channels follow from these earlier results upon a series of replacements [T7] . 
This means that the twist-3 light-front distribution amplitudes are assumed to be determined by 
two-particle configurations only. The amplitude for —>■ foK^ decay is given in Ref. [T7j, so that 



we identify, after Eq. (10) 



An 



+ 



G p 
G p 
G p 



<-^X<-oKo-^f«8) 



foK 



fKF^f^\Ml,){Ml-M 



:i2) 



foK 



fKF^f\Ml,){Ml-M 



+ 



A' 



+ 



+ 



G p 
71 
G p 
71 
G p 
71 

Gp^ 
V2 

G p 
71 



Jb , 



(13) 



^6 - o«S 



Kfo 



r 



«6 - 



Kfo 



^ / Kfo 



f/;^Fo^^(Mj,)(M|-M^o) 



(14) 



(15) 



The M1M2 subscripts mean that the quantities in brackets are to be interpreted as af(MiM2), in 
which the Mi meson contains the spectator quark from the S-meson, and bi{MiM2) , in which the 
Ml meson carries an anti-quark from the weak vertex and the M2 meson contains a quark from 
the weak vertex. The QCD factorization framework is rigorous in leading order in 0{AQCD/fnb), 
but it also includes estimates of 1 /mf,-suppressed corrections. We note that f^" = 2Mf^/mb{iJi) 
and = 2M|,o/("^fe(At)("^s(At) + fndifJ'))) but that these chirally-enhanced terms are counted as 
terms of leading power in the l/rrif, expansion in QCD factorization. The quark masses are running 
masses defined in the MS scheme. Certain l/rrib corrections can suffer endpoint divergences, and 
their estimate is uncertain. The 6j(MiM2) terms, for example, reflect annihilation contributions, 
which are a class of l/m^-suppressed corrections to the decay rate. Moreover, the a^{MiM2) terms 
contain 0{as) corrections to the operator matrix elements, which can, in turn, contain l/m^- 
suppressed corrections. The 0{as) corrections encode non-perturbative input through integrals 
over the light-cone distribution amplitudes. The coefficients af(MiM2), where p = u,c, can be 
expressed as 



af (M1M2) = Qifi) + 



Nr. 



+ 



V,{M2) + —H,{MiM2) +Pf(M2), (16) 
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where the upper signs apply when i is odd, the lower signs apply when i is even, and Cp = 
(iV2-l)/2iVc with = 3. The quantities ^^.(Ms), //^(MiMa), and Pi(M2) reflect vertex corrections, 
hard spectator interactions, and penguin corrections, respectively. The function Hi{MiM2) contains 
a endpoint divergence in its power suppressed terms — we refer to Refs. [291 EO] for all omitted 
details. For work towards a theory of the power corrections, we refer the reader to the developments 
in Refs. pT| [55] . The coefficients 6j relevant to our calculation can be expressed as 



53(MiM2) 



Cs A\ + C,{Ai + 4) + C, Ai 
'c,A', + C,{A'^ + A{) + N,Cs4 



(17) 



where in the annihilation amplitudes A][^ = y4^-^(MiM2) the superscripts i and / refer to gluon 
emission from initial- and final-state quarks, respectively. For the calculation of A^^^ {M1M2) we 
include the corrections coming from a^, the second Gegenbauer moment of the kaon, for consistency 
with our analysis of the /o(980), and we provide explicit expressions in App. |Aj 

The expressions for B f^K decay contain non-perturbative hadronic input through the meson 
decay constants and B M form factors; these quantities are sensitive to the assumed quark 
structure of the hadrons. Assuming the /o(980) resonance can be written as a qq state, we write 

|/o(980)) = cos^lss) + sin^lnn) , (18) 

where \nn) = {\uu) + \dd)) /\/2. The empirical observation of both r(J/-?/; /quj) andT{J/^ — > /o0) 
suggests that the /o(980) has both strange and non-strange components, so that 6 is non-zero [16j . 
We define the B ^ fo form factor as 



{foip)\di,i,b\Bip)) 



Ml 



JO „ TP^JOl 



r J q 

where P^ = {p + p')^ and q^i = {p ~ p')^i, and the /o decay constant as 

{fo\qq\0) = MfJl 
for q G (n, s). Defining |/q) = \qq), we have 

1 



{f^\ss\Q) = Mf,f% and {f^\uu\0) 



V2 



fo J fo ' 



(19) 



(20) 



(21) 



so that ff^ = ff^ cose and ff^ = /j? sin^. Similarly we define Fq^° = smOF^^^ where F^^° 
describes the form factor to the \dd) piece of the /o(980) final state. 

In order to calculate the B fo form factor we assume the fo to be a qq state and use the 
constituent quark model (CQM) of Refs. [IHl HH 1121 IS], which combines heavy quark effective 
theory with chiral symmetry in the light quark sector. In Refs. |10l HI], Gatto et al. study the 
D — > (TTT — *• Stt and Dg —>■ /ovr amplitude using the CQM model for the D —> a and Ds — > fo form 
factors and find good agreement with E791 data. This model describes interactions in terms of 
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effective vertices between a ligfit quark, a lieavy quark, and a lieavy meson. Tlie model depends 
on botli its UV and IR cutoffs. Tlie UV cutoff A is set by tlie spontaneous cliiral symmetry 
breaking scale A^, which is of (9(1 GeV). This model does not include confinement, so that one 
has to introduce an IR cutoff /2. The constituent quark mass m, is determined by solving the 
Nambu Jona-Lasinio (NJL) gap equation, which, in turn, depends on the UV and IR cutoffs. We 
choose A = 1.25 GeV as described in Ref. |13]. Thus for fixed A, as the IR cutoff varies, m varies 
accordingly, as explicitly illustrated in Ref. For the default value of the B fo form factor 
we use m = 0.3 GeV, Ji = 0.3 GeV, and Ah = 0.4 GeV. The parameter Ah = Mh — Mq, where 
Mh is the mass of the heavy meson and Mq is the mass of the constituent heavy quark. Explicit 
expressions for the polar and direct contributions to the form factor are given in Refs. [101 E] • To 
assess the range of the B ^ fo form factor we vary fi in the range [0.25, 0.35], as per the associated 
variation in m given in Ref. [H], and vary Ah in the range [0.3, 0.5] as given in Refs. |10l UH Il2| B3]. 
We find that the variation in F^-^o is mainly driven by the variation in /2. However, the error range 
reported in Refs. |Sni HI] was apparently determined by varying Ah alone [45]. The default value 
of the B ^ fo form factor, as well as its uncertainty, are given in Sec. Ill A[ 



The uncertainty in the calculation of the decay amplitude comes from both the statistical and 
systematic errors present in the QCD factorization approach. The uncertainties associated with 
different input parameters such as the scalar meson decay constants, the form factors, the quark 
masses, and the Gegenbauer moments of light-cone distribution amplitudes are of first kind, whereas 
the uncertainties associated with the 1/mb-suppressed corrections are of the second kind. The the- 
oretical uncertainties in the latter case come from the hard spectator and the weak annihilation 
contributions which contain endpoint divergences. The endpoint divergences Xh in the hard spec- 
tator and Xa in the annihilation terms are parameterized as ||29j 



Xh,a = In (1 + PH,Ae"^"-^) , (22) 

where we assume ph,a < 1 and = 0.5 GeV. We note 4>h,a are unknown, strong-interaction 
phases. 

We wish to determine the impact of the various theoretical uncertainties in the B f^K decay 
amplitude on the value of ASf^^Xs ^ quantitative way. To realize this, we begin by defining a 
"default model." This consists of using the central values of the inputs given in Ref. [T71 129], as 
well as setting p^i = p^^ = in the parameterization of the endpoint divergences. With this in place 
we thus determine the default values of the A"'^ amplitudes. To gauge the size of the uncertainties, 
we perform a random scan of the allowed theoretical parameter space. That is, we include the 
uncertainties coming from all the input parameters. For the parameter scan, we choose the range 
of the input parameters by taking either la, 2cr, or 3cr deviations from the central values. Once we 
have the maximal and minimal values of the inputs in the chosen ranges, we draw random values 
of the input parameters within that range. Similarly, to include the uncertainties coming from Xa 
and Xh we vary ph,a from to 1 and (pH,A from to 27i. That is, we draw random values of ph,a 
in the range to 1 and 4>h,a in the range to 27r. 
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III. RESULTS AND DISCUSSION 



A. Inputs 

For definiteness, we summarize the parameter choices of our default model. We employ a renor- 
malization scale of 0{mh/2); namely, we choose /i = 2.1 GeV, so that the value of the strong cou- 
pling constant in NLO at this scale is 0^(2.1 GeV) = 0.286. To realize this, we work with Nf = 5 
throughout, after Ref. [29j, and choose A*^^) = 0.225 GeV, which corresponds to as{Mz) = 0.118. 
As a check of the accuracy of NLO precision in this context, we compare our results with those 
computed at /i = 4.2 GeV, for which we note as{mi,) = 0.224. For the electromagnetic coupling 
constant, we use = 129 and neglect the evolution of a as in Ref. p9]. To realize the Wilson 
coefficients, we follow the procedures of Ref. [29], using mtirrit) = 167 GeV, Mw = 80.4 GeV, and 
siia^ 9w = 0.23 j^Sj ^ and verify the leading-order (LO) and NLO Wilson coefficients they report 
through exphcit computation. We also use the values = p// = to fix the default values of the 
hard spectator and annihilation terms. We fix the value of the b quark mass at the scale p = 2.1 GeV 
to be m6(2.1GeV) = 4.88 ± 0.08 GeV, which corresponds to mb{mb) = 4.2 ± 0.07 GeV |16] us- 
ing the two-loop expression for the running quark mass [M]. The error in mf,(2.1GeV) is cal- 
culated using the maximum and the minimum value of rUb at the p = rrib scale. Similarly the 
value of ms{mb) = 0.077 ± 0.017 GeV corresponds to m,(2.1GeV) = 0.090 ± 0.020 GeV, after 
Ref. [30]. For the meson masses we use Mb = 5.2795(5) GeV, M^o = 0.497648(22) GeV, and 
= 0.980(10) GeV, as given in Ref. [IS], where the uncertainty in the last digits is indicated by 
the number in parentheses. Similarly, the mean life of B meson is tb = 1.530(9) x 10~^^ s ^6], 
and Gf = 1.16637(1) x 10"^ GeV'^ [16]. The errors in these empirical quantities are unimportant 
for our purposes, so that we neglect them in our random scan. The CKM matrix elements, taken 
from [17] are, 

IKJ = 0.00357+°-°°°|^, IKJ = 0.22653+°-°°°^^ 



ub\ — "•'J'J"^"-' ' —0.00017' {rusl — "•■^^'J"'-'— 0.00077 

M — n n/inc;+o-oo32 it/ \ _ n q7qi«+o.oooi8 

cb\ — U.U^UO__Q QQ29, \Vcs\ — U.y(0±U_QQQQ]^g 



\Xu/Xc\ = 0.0205 ± 0.0019, |Ac| = 0.0394 ± 0.0031 

7 = 76.8^^?:^ /3 = 21.51};°, (23) 

where 7 and (3 are reported in degrees. We use the uncertainties associated with the CKM matrix 
elements to calculate the uncertainties in |Ac| and |A„/Ac|, assuming the errors are uncorrelated for 
simplicity. We note the angle 7 is the phase associated with Xu/K, namely A^/Ac = lA^j/Ad e~^^, 
and that f3 is the unitarity-triangle angle we defined in Sec. [T] In our scan over parameter space we 
vary our inputs within 2a and 3a of their default values as well. In the case of the CKM parameters 
we use the ranges as reported in Ref. [47] for |Vifc|,7, and /? for ±2a and ±3a as appropriate; it is 
worth noting, in particular, that 7 when ranged over a 3a variation is never negative. 



Repeating our calculation with rut ~ 172.7 GeV, as per Ref. jJS], and following the procedures described in text 
yield ASf^Ks = 0.0269, Cf^Ks = -0.00557, and Br{ foKs) = 13.6 x 10"^ in place of the values reported in Eq. (|31 1. 
We can safely neglect this update in our parameter scans. 
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The other inputs are taken from Refs. [T71 [291 SH] and are given as follows. All the scale- 
dependent quantities in the scalar sector are evaluated at /i = 2.1 GeV as per Ref. [17]. The ratio 
of charm quark mass to the b quark mass is taken from Ref. |1H]; the resulting value of the charm 
quark mass encompasses the value recommended in Ref. [16]. The value of the mixing angle is taken 
from Ref. [17]. All the other input parameters are taken from Ref. [29]. 

5i = -0.54 ±0.06, S3 = 0.01 ±0.04, 

af = 0.3±0.3, af = 0.1±0.3, 

Xb = 0.35 ± 0. 15 GeV , /b = 0.20 ± 0.03 GeV , 

Fq^^ = 0.35 ± 0.03 , /k = 0.16 GeV, 

ff, = 0.460 ± 0.025 GeV , 9 = 152.5 ± 12.5° , 

mc/mb = 0.27 ±0.06, mg/m, = 0.0413 . (24) 

We note that rriq = {rriu + md)/2 and that we neglect the error in mq/uis. This amounts to 
neglecting the error in the light-quark mass m^, since we include the error in the strange quark 
mass. This we may safely do as we employ mg/m^ in the evaluation of exclusively; we note 
that mass differences of 0{md — m„) are tantamount to the inclusion of isospin-breaking effects, 
which are numerically unimportant to us here. Explicit expressions for the leading-twist, light-cone 



distribution amplitudes of the light mesons are given in Eq. (Al), where we note Bi and B^, are 
the first and third Gegenbauer moments of the /o(980) and af and are the first and second 
Gegenbauer moments of the kaon. Although all the Gegenbauer moments are scale-dependent in 
principle, the inclusion of such variations is beyond the accuracy of the NLO treatment we effect 
here, so that we neglect such refinements. The first inverse moment of the B meson light-cone 
distribution amplitude. As, is defined by dx^B{,x)/x = Mb/Xb, where is one of two 

twist-2 light-cone distribution amplitudes of the B meson [291 [30]. We term /b, Fq^, /k, and ff^ 
the B meson decay constant, the B ^ K form factor, the K meson decay constant, and the scalar 
decay constant, respectively. The given value of is specific to the assumed qq structure of /o; 



we recall 6 is the mixing angle of Eq. (18). For our inputs, we note that the ratios for the K 
meson and /o meson, are 

= 1.08 , f^" = 0.402 . (25) 

Our default value of the B fo form factor is 

Fq^° = 0.284 , (26) 

where we use the CQM model and the parameters m = 0.3 GeV, jl = 0.3 GeV, and Ajj = 0.4 GeV. 
Cheng et al. [T^ assert that the F^-^o form factor should be comparable in size to the B —>■ n form 
factor, and we find this to be consistent with our own numerical estimate. We calculate the error 
in the S — /o form factor by varying fx in the range [0.25, 0.35] and Ah in the range [0.3, 0.5]. We 

Ft f ^ 

find that the variation in Fq ° with respect to A// is small compared to the variation with respect 
to ft. That is, varying Ah over our chosen range makes Fq ° range over the values [0.28,0.29]. 
Varying fi over its chosen range as well, we find that Fq ° ranges from [0.23, 0.46]. For the random 



10 



scan we use the range [0.23, 0.46] for our la scan and triple its range, to yield [0.0048, 0.69] for our 
3a scan. 



The coefficients a^{foK), a^glKfo) and 6j are calculated at the scale /i = 2.1 GeV using Eqs. (16) 



and (17), as well as the formulae in Ref. [29], to yield 

< = -0.0261 - i0.0208, < = -0.0340 - zO.OllO, 

a;f = -0.0581 - i0.0185, = -0.0641 - i0.00842, 

= (76.6 - iO.434) x 10'^ = (76.5 - zO.259) x 10-^ 

= (-172 + zl30) X 10^^ a^o = (-172 + il30) x 10-^ 
als{Kfo)=al,ifoK), 

hifo Ks) = -0.0506, 63(^5 /o) = 0.0264, 

b3,Ew{fo Ks) = -0.00133, bs,Ew{Ks fo) = -0.000768, (27) 

Our formulae for 63 and b^EW differ slightly from those given in Ref. [HI |29], as we detail in App. [A| 
though the numerical differences are negligible. Now that we have all the input parameters, we can 



calculate the default value of the A^'^ amplitudes using Eq. (12). We find 

Al = (-8.12 + zO.438) x 10"^ GeV, A"^ = (-7.67 - i0.821) x 10"^ GeV, 

= (-84.7-i24.2) X lO^^GeV, = (-92.5 - ill.O) x 10"^ GeV, (28) 

If we ignore the annihilation terms from the B foKs decay amplitude, we find 

Al = (-12.5 + zO.438) x 10"^ GeV, A^ = (-12.1 - i0.821) x 10~^ GeV, 

Al = (-76.4 - i24.2) x 10"^ GeV, A^ = (-84.3 - ill.O) x 10~^ GeV, (29) 

If we ignore all the l/rrib power suppressed terms in the B — > foKs decay amplitude, we find 

Al = (-12.1 + ^0.438) X 10-^ GeV, A^ = (-11.7 - i0.821) x 10"^ GeV, 

Al = (-76.4 - i24.2) x 10~^ GeV, A^ = (-84.3 - ill.O) x 10~^ GeV, (30) 

For our default parameter set, the amplitudes mediated by the ss component of the /o(980) are 
numerically dominant, and the various contributions coming from the power-suppressed terms are 
small corrections to the leading power contributions. However, these power corrections suffer end- 
point divergences, and their estimate is uncertain; they could well affect the value of ASf^Ks- 



shall consider this possibility in detail in Sec. IIIB We note that the dominance of the penguin 



contributions associated with the strange quark component of the fo emerges because the 04 and 
Og interfere destructively exclusively in the decay to the nonstrange component of the /o, as we can 



see from Eq. (12). Moreover, we note that the contribution coming from the annihilation ampli- 
tude in this case is smaller than the contribution coming from the penguin amplitude, so that the 
S — > /o Ks decay amplitude is dominated by the ss component of the /q. 



B. Two-body Decay Results 

Using the default values of the A"'^ amplitudes we can easily calculate df^^Ks ^-nd thus ASf^Ks 
and CfgKs using Eq. ^ and Eq. ([7]). We compute the values of AS'/^Xs, ^/o^s' branching 
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ratio (Br) for the two-body decay as discussed in Sec. |TTj For our default parameter set we find 

ASf.Ks = 0.0269, Cf.Ks = -0.00561, Bi{foKs) = 13.4 x 10"^ . (31) 

We did the same calculation ignoring the annihilation terms, for which we find 

ASf.Ks = 0.0268, Cf.Ks = -0.00584, Br{foKs) = 12.4 x 10"^ . (32) 

Alternatively, if we neglect all l/rrib suppressed terms we find 

ASf.Ks = 0.0268, Cf.Ks = -0.00587, BiifoKs) = 12.3 x 10"^ . (33) 

We have done the same calculation at the fi = rrib scale, and our results are very similar to those 



given in Eq. Q, Eq. (Q, and Eq. Q. That is, we find the values of ASf^Ks are 0.0269, 0.0269, 
and 0.0269, respectively, whereas the values of Cf^Ks are —0.00603, —0.00621, and —0.00623, 
respectively. However the value of the branching ratio does change, giving 10.3 x 10~^, 9.69 x 10~^, 
and 9.63 x 10^^, respectively. It is evident from our numerical estimates that ASj^Ks and Cf^Ks 
and the two body branching ratio receive little contribution from the annihilation terms or, indeed, 
from any of the l/nib suppressed terms. That is what we expect since the B —> foKs decay 
amplitude is driven by the s-quark component of the /o; the leading power contributions to the ^ 
amplitudes are much larger than the contributions coming from the power suppressed terms, and 
ASfgKs and Cj^Ks depend only on the ratio of the A^'* amphtudes df^Xs- However, the branching 
ratio is controlled by the square of the sum of the amplitudes. Thus any change in the size of the 
amplitudes themselves will definitely change the branching ratio; this explains why our branching 
ratios are more sensitive to the value of fi. Our default value of the branching ratio is not consistent 
with the experimental value of (5.8 ± 0.8) x 10~^ obtained by BaBar and Belle [19], ED], so that 
the default set of input parameters we employ can not explain the experimental data. This does not 
falsify our approach per se, as it is reasonable to think that we can compute the ratio of amplitudes 
df^Ks niore accurately than the A^'^ amplitudes alone. We also find that we are able to confront 
the empirical branching ratio successfully in our scan of the theoretical parameter space. 

The numerical values of ASf^Xg and Cj^Ks for B foKs decay are calculated by Ref. [E] as 
well, for which they find 

ASf.Ks = 0.023, Cf.Ks = -0.008 . (34) 

Our default values of ASj^^Ks and Cf^Ks for B — > foKs decay are similar to those in Ref. [T7], 
but they are not exactly same. We note that we also differ in our computed values of a^{foK) and 
a^g^Kfo), though this emerges, at least in part, because they employ a value of q:s(/x = 2.1 GeV) 
appropriate to Nj = 4. 

Our default values of ASf^^Ks and Cf^^Ks are small and that is what we expect from the SM 
point of view. In the QCD factorization approach the uncertainties in the calculation of the decay 
amplitude come from the input parameters as well as from the l/rrib suppressed terms. We would like 
to establish the range of ASf^Xs and Cf^Xs possible from SM physics using the effect of the various 
uncertainties in the B — >■ foKs decay amplitude. Thus it is necessary to include the corrections 
arising from the uncertainties associated with all the input parameters, including those in the scalar 
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meson decay constants, the form factors, the quark masses, the CKM elements, and the Gegenbauer 
moments of the hght-cone distribution amphtudes. Since the quark structure of the /o is not well 
known, we are particularly interested in whether the uncertainties in the B —>■ fo form factor and 
the /o scalar decay constant impact the value of ASf^Xs ^ significant way. Additional systematic 
uncertainties in our calculation come from the power corrections which contain endpoint divergences; 
our estimate of these is uncertain. As we have noted, too, the branching ratio we compute from 
our default parameter set does not confront the experimental branching ratio successfully, though 
this is not necessary to describe ASj^Ks ^^1^- 

To see the effect of the above mentioned uncertainties on the various observables, we perform a 
random scan of the allowed theoretical parameter space. The central values of the inputs and the 
uncertainties associated with them are in Eq. (23), Eq. (24), and Eq. (26). For the parameter scan, 
we choose the range of these input parameters by taking la and 3cr variations from the central 
values. Once we have the maximal and minimal values of the inputs in the range mentioned above, 
we draw random values of the input parameters in the chosen range. To include the uncertainties 
coming from the hard spectator and the weak annihilation terms we vary Ph,a from to 1 and 4>h,a 
from to 27r. That is, we draw random values of Ph,a in the range to 1 and (f)H,A in the range to 
2n. We use the public domain random number generator "rannyu" throughout; note Ref. [51] for 
a discussion of suitable inputs. Each chosen parameter set corresponds to a theoretical model, and 
we have plotted the various combination of ASf^Ks ^/o^s ^^0' models for a particular 
seed, just for illustration, for parameter ranges fixed at la and 3a from their central values in Fig. [Tj 
The box-like shapes of the resulting regions suggest that there is little correlation between ASf^Ks 
and CfQKs- That means {df^xsl is small, and ASf^Ks ^foKs mainly driven by Re^df^Ks) 
and Im((i/Qi^g), respectively, as we can see from Eqs. ^ and ([T]). 

We did a random scan for two different seeds and for each seed we use 100, 000 and 500, 000 
different theoretical models. We determine the range in AS/^Ks which results, or indeed of any 
observable, by evaluating the extremal values which emerge from the parameter scan. This makes 
the results sensitive, in principle, to the detailed manner in which the scan is effected. The resulting 
range in ASf^Xs varies little for different values of the seed and for the different numbers of theo- 
retical models if we choose the input parameters within la of the central values. For example, the 
range of ASf^Ks ^ scan is found to be [0.017, 0.034] and [0.016, 0.035] for two different seeds 
from a sample of 500,000 parameter sets. However, the range of ASf^^Ks ^'^^^ change significantly 
with seed if we choose the input parameters to range within either 2a or 3a of the central values. 
For example, the range of ASj^Ks scan over a 3a variation is found to be [—0.41, 0.29] and 

[—0.34, 0.20] for two different seeds from a sample of 500, 000 parameter sets. It is evident from our 
scan that though the range of ASj^Ks small for la variations, it is large for 3a variations. No 
color-suppressed tree amplitude contributes to the B foKs decay process, but we have large ex- 
cursions in ASf^Ks small regions of the parameter space nevertheless. In these special regions, 
the amplitude becomes small and drives the large values of ASf^^Ks- However, such small values 
of always give theoretical branching ratios which are much too small. It thus becomes crucial to 
apply the experimental branching ratio constraint to determine the allowed range in ASf^Xs within 
the SM. We recall that the average branching ratio for B foKs decay process measured by Belle 
and BaBar [IHl ED] is (5.8 ± 0.8) x 10~^. We impose the branching ratio constraint in such a way 
that we ignore those theoretical models which are not compatible within la of the experimental 
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branching ratio for la parameter scans and within 3a of the experimental branching ratio for the 
3a parameter scans, respectively. After we impose the empirical branching ratio constraint, the 
seed-averaged range of ASf^Xs scan of models spanning la and 3a variations are found to 

be [0.018,0.033] and [—0.019,0.064], respectively, for the 500,000 point simulation. In this context 
"seed-averaged" means that we report the extremal values of ASf^Ks which result from scans using 
two different seeds. It is interesting to note that the range in ASf^Xs oiice the branching ratio con- 
straint is applied is almost identical to what was found without the branching ratio constraint in 
the la scan. This underscores our point that the branching ratio is largely independent of the value 
of ASf^Ks- ^foKs well can range over a wide array of values giving [—0.47,0.51] 

and [—0.55, 0.49] for two different seeds from a sample of 500, 000 parameter sets within 3a of the 
default values. However, after applying the branching ratio constraint, such large excursions disap- 
pear, giving the seed averaged range [—0.045, 0.051] for 3a and [—0.013, —0.0024] for la. We plot 
the various combinations of ASf^Ks ^foKs with the branching ratio constraint for the seed 
used in Fig. [T| just for illustration, for a random sample of 500, 000 models in Fig. |2] After applying 
the branching ratio constraint, we find that the range in ASf^Xs ^foKs depends on neither 

the number of parameters nor the particular seed used in the random scan. To test the accuracy 
of our NLO analysis, we also change the renormalization scale /i from mi,/ 2 to m^, and we find 
very little variation in ASf^Ks ^^"^ ^/o^s oiice the empirical branching ratio constraint is imposed. 
That is, the range of ASf^^Xs ^/o^s ^^^^ '^^^^ [0.017, 0.032] and [—0.012, —0.0032] for la, 
respectively, whereas they are [—0.015, 0.061] and [—0.044, 0.033] for 3a. The weak /i dependence 
we observe follows as the observables involve ratios of computed amplitudes in each case. 
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FIG. 1: Range in ASf^Xs ^^^^ ^foKs froni a scan of 500,000 theoretical models. The lighter interior 
region corresponds to a scan of the parameter space at la, where the darker, larger region corresponds to 
a scan at 3fT. Some points in the 3a scan fall outside the window chosen for the illustration. 

To study the likelihood of various values of ASj^Ks theory space, we have plotted a histogram 
of ASjqKs from a sample of 100,000 and 500,000 theoretical models, after averaging over seeds, in 
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FIG. 2: Range in ASf^Xs ^foKs from a scan of 500,000 theoretical models with the empirical 
branching ratio constraint imposed. The darker points represent the possible range of ASj^Ks ^/o^s 
within 3o" whereas the lighter points represent the range for a la variation. 

Fig. [3j To plot the histogram we first set a bin size and then determine the number of models 
which fall within each bin. Since we have different numbers of theoretical models in our random 
scans, to plot all of them together we divide the number of models falling in each bin by the total 
number of theoretical models used and choose that quantity as the ordinate of each histogram. If 
no points are shown, both here and in later figures, that means that no models whatsoever occupy 
that region of theory space. For scans employing parameters within a la variation the range of 
ASfgKs varies little if we change the seed or the number of models used in the scan. However, it 
does vary notably if the parameters chosen range within 2a or 3a variations of their central values. 
In Fig. |4| we plot a histogram of ASf^Xs: after averaging over seeds, from a sample of 100, 000 and 
500, 000 models taking the values of the input parameters within 3a of their central values. It can 
be seen from the histogram that although the actual range of ASf^^Ks varies for different seeds and 
for different number of parameter sets, the shape of the histogram in ASfgKs does not seem to vary. 
We also plot the histograms of ASf^Xs after applying the branching ratio constraint, as shown in 
Fig. [5] and Fig. [6} and note that the shape of the histogram changes little. These results are also 
seed-averaged. In these figures the total number of points refer to the total number of models which 
satisfy the branching ratio constraint. In particular, the scale of the ordinate is not same in Fig. |5] 
and Fig. |6] because the total number of points which satisfy the branching ratio constraint at la is 
small compared to those which satisfy the branching ratio constraint at 3a. Note that the integral 
of these model space results over ASj^Ks should yield unity in each case. 

We have also studied the impact of the hadronic uncertainties alone on ASf^Xs ^'^^ ^ ~^ fo^s 
decay process. To do this we set the CKM parameters to their default values and vary all the other 
input parameters. We include the variations in the power corrections as well. Similarly to see the 
impact of the CKM parameters, we set all the other hadronic inputs to their default values and vary 
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FIG. 3: Histogram for ASf^^x^ random scan employing parameters which range over la of their 

central values. We note + denotes the scan with 100,000 models and x denotes the scan with 500,000 
models. The range of ASj^Ks same for each case. 
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FIG. 4: Histogram for ASf^Ks ^ random scan employing parameters which range over 3a of their 
central values. We note + denotes the scan with 100,000 models and x denotes the scan with 500,000 
models. Some points in the scan fall outside the window chosen for the illustration. 
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FIG. 5: Histogram for ASj^Ks^ after that of Fig. [3j once the empirical branching ratio constraint is 
imposed. Here the total number of points refer to the total number of models which survive the imposed 
branching ratio constraint. 
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FIG. 6: Histogram for ASj^Ks^ after that of Fig. |4j once the empirical branching ratio constraint is 
imposed. Here the total number of points refer to the total number of models which survive the imposed 
branching ratio constraint. 
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only the CKM parameters |Ac|, |A„/Ac|, 7, and /3. We include the uncertainties associated with Vij 
in an uncorrelated way to find the uncertainties in |Ac| and |A„/Ac|. We perform a random scan over 



these CKM parameters within la and 3a of their central values as given in Eq. (23) and Ref. |47j . 



The associated histograms of ASf^Xs scan are shown in Fig. [7] and Fig. [8] These results 

are also seed averaged. As we see from the figures, the impact of the hadronic uncertainties are 
very small compared to the CKM uncertainties in this case. The actual shape and range of ASj^^Ks 
is driven mainly by the CKM parameters. Since \dfgKs\ is small, the second term in the numerator 
of Eq. ^ is small compared to the first term — ASj^Ks mainly driven by the first term in the 
numerator. Both terms do tend to zero, however, as 7 becomes small. Thus we expect to get 
very small values of ASf^Xs sufficiently small values of 7. As we can see from the histograms, 
variations in 7 and, indeed, the CKM parameters impact the likely values of ASj^Ks- Since 7 is 
never negative and the \dfgKs\^ term is negligible, AS* can only be negative if Redf^^Ks < 0, which 
is an exceptionally rare occurrence. 
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FIG. 7: Histogram for ASf^Xs varying only the CKM parameters over a la variation. 



When we repeat this analysis for a scan in which the parameters are allowed to range over 3a, 
as shown in Fig. |9] and Fig. 10 , we find that the hadronic uncertainties dominate over the CKM 



uncertainties. This emerges, despite the detailed shapes in Fig. |9] and Fig. [T0| because we use the 
extremal values found in the simulation to define the range. The range of ASj^Ks thus found 
to be [—0.34, 0.25] for the hadronic uncertainties and [0.00076, 0.039] for the CKM uncertainties. 
However, after applying the branching ratio constraint the ASf^Xs range is found to be small, 
giving [—0.018, 0.048] for the hadronic uncertainties and [0.00080, 0.038] for the CKM uncertainties. 
In the 3a scan the range of ASf^^Xs is large, so that \dfgKs\ is no longer small, for certain values 
of the parameter set. It is the small value of the P'^ amplitude which is responsible for these large 
values of Idf^Xsl- For A* = "^b, the ranges of ASf^^Xs quite similar giving [—0.0073,0.043] for 
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FIG. 8: Histogram for ASf^Xs varying only the hadronic parameters over a la variation. 

the hadronic uncertainties and [0.00075,0.039] for the CKM uncertainties once the branching ratio 
constraint is imposed. 
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FIG. 9: Histogram for ASf^^Ks varying only the CKM parameters over a 3a variation. 
We have also studied the impact of the individual hadronic parameters on the value of ASf^^Xs 
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FIG. 10: Histogram for ASf^^s varying only the hadronic parameters over a 3a variation. More points 
exist beyond the window chosen for the illustration. 

for the B foKs decay process. To do this we set all the other input parameters to their default 
values except the one which we want to study. Note that we include the variations in the power 
corrections as well throughout. We also use the same seed used to generate Figs. [T]and [2j We find 
that for any hadronic input with values in the la range, the range of ASj^Ks very small. For the 
input parameters with values within the 3a range, the largest excursions in ASf^Xs come from the 
first inverse moment of the B meson distribution amplitude A^, which enters the B —>■ JqKs decay 
amplitude through the hard spectator interaction terms. However, after applying the branching ratio 
constraint the ASf^Ks range is found to be small giving [0.026, 0.034] if we employ a renormalization 
scale of /i = 'mi,/2 and giving [0.026,0.034] for /i = m^. The negative values of ASfgXs found in 
the scan with inputs which range up to 3a of their central values are driven by the uncertainties 
associated with A^, Xa, and Xh- To see the impact of these parameters alone on ASf^Xs^ 
the other input parameters to their default values and vary A^, Xa, and Xh within 3a from their 
default values. The range of ASj^^Ks quite large giving [—0.25, 0.19]. However, we find ASj^^Ks ^'^ 
be small and positive once we impose the branching ratio constraint, giving the range [0.024, 0.036] 
and [0.025,0.036] for the two different renormalization scales we employ, /x = mb/2 and /i = mb, 
respectively. Fixing the values of A^, Xa, and Xh to their default values, we perform a random 
scan over the other input parameters within 3a and find the range of ASj^Ks t)e [—0.0015, 0.047]. 
Once we impose the branching ratio constraint we find the range of ASj^Ks ^'^ ^® [0.00054, 0.046] 
and [0.00037,0.046] for = 'mb/2 and ^ = rrib, respectively. It is interesting to note that the range 
of ASjqKs fi^*^ t)y fixing the values of A^, Xa, and Xh to their default values, and imposed the 
branching ratio constraints, captures most of the range we find once we perform the scan with all 
the parameters within 3a. However, we note that the extension of those scan results into negative 
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values come from the variations associated with A^, Xa-, and Xh as welL 

In principle, ASjp Ks can have large excursions either by having a small P'^ amplitude or a large 
amplitude. The first possibility can be controlled by demanding that the theoretical model confront 
the empirical branching ratio successfully. As we have already noted, such small amplitudes give 
theoretical branching ratios which are much too small, so that they can be excluded by imposing 
the experimental branching ratio as a constraint. We now wish to consider whether it is possible 
to have a large P^ amplitude without enhancing P'^ as well. If it were possible, one could have a 
large value of ASj^^Ks without necessarily violating the branching ratio constraint. In the case of 
the default parameter set, the a\, Og, and a^Q — a\ is so small as to play no role — terms interfere 
destructively giving small A^^ amplitudes. In this case, the hard spectator terms, which enter 
through the a\ and a^Q terms, are similar in magnitude to the vertex terms in Og. However, once 
we do a random scan over the complete theoretical parameter space, the hard spectator terms play 
a dominant role in enhancing the a\ and a^o terms, which, in turn, enhance the A^^ amplitudes. 
The real and imaginary parts of a\ — O.Sa^Q, noting Eq. 12, e.g., can both be very large, so that 
the cancellation of these terms with Og is no longer possible in these regions of parameter space. 
However, these effects act to enhance both and P'^ at the same time — it is not possible to 
have a large P" amplitude without having a large P^ amplitude as well. Since both A^ and A^ 
amplitudes are large, ^Sf^Ks is not enhanced. Thus large excursions in the hard scattering terms 
can not produce large excursions in the value of AS'/qJ^^. 

Our primary motivation in conducting our current analysis is to determine whether it is possible 
to control the range of ASf^^Ks irrespective of the structure of the /o resonance. In Sec. II we noted 
that the structure of the /o impacts both the B —>■ Jq form factor and /o decay constant. Thus 
the impact of these parameters alone on ASj^^Ks gives us insight as to whether we can constrain 
the range of ASf^Xs small values no matter the structure of the /q. We perform a random 
scan varying only these two parameters within la and 3a from their central values, both with 
and without the branching ratio constraint, for the particular seed used to generate Figs. [T]and[2j 
Without the branching ratio constraint, the range of ASj^Ks is [0.026, 0.028] for la and [0.025, 0.028] 
for 3a, respectively. Once we impose the branching ratio within 3a of the experimental value as 
a constraint, we find the range to be [0.025, 0.028] for the scan in which the B —>■ fo form factor 
and /o decay constant are varied within 3a of the central value. For = mi,, the range of ASf^Ks 
in this case is also very small, giving [0.026,0.028]. It is evident from our analysis that the impact 
of the B ^ fo form factor and fo decay constant on ASf^Xs is small irrespective of the branching 
ratio constraint. Apparently the strange quark content of the fo is so important in determining 
ASfgKs ^^^^ particular value of the ratio of the B ^ fo form factor and /o decay constant is 
quite unimportant. 



IV. CONCLUSIONS 



In this paper, we have analyzed the two-body decay of the B meson to one scalar meson and 
one pseudo-scalar meson, namely to the fo(980) Ks final state. The empirical study of this mode 
measures sin(2/3) through a decay process controlled by the b —>■ s penquin amplitude. This is 
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important for several reasons. First, the study of such modes is complementary to the study of 
the 6 — i> s transition in Bg mixing. The hint that the empirical weak phase determined from Bg 
mixing is larger than the SM prediction [52], if borne out, would require new sources of CP violation 
beyond that contained in the CKM paradigm. This makes the confirmation of, or constraints on, 
such effects in the b s penguin modes crucial. Generally, we would like to assess the possible 
deviation of Sf{t), the time-dependent CP asymmetry, from SM physics in a robust way, so that we 
can determine what new physics, if any, exists in these decay modes. To that end, the /o(980)i^5 
final state itself is of intrinsic interest. The time-dependent CP asymmetry in this case has no color- 
suppressed tree contribution; it is intrinsically subject to less theoretical uncertainty than modes in 
which such contributions are present. Moreover, its empirical study can be complementary to that 
of another theoretically clean process: B (pKg, as both modes occupy the B KKKg Dalitz 
plot. Indeed, B —>■ f^Kg and B <f)Ks, assuming the (p to be ideally mixed, are the only known 
b s penguin modes which have no color-suppressed tree contributions. 

We have investigated the size of ASj^Ks^ i-^-' deviation of Sf(t) from sin(2/3), in the B 
foKs decay process in the SM using the QCD factorization approach, assuming the /o to be a 
qq state. We employ a parameter scan to probe a broad range of possible theoretical models, 
exploring variations in the inputs at the 3a level and ill-known 0{Aqcd/Mb) corrections with 
100% uncertainty. The B — > foKs decay mode has been studied by other authors within the 
QCD factorization approach [151 UHl E] • Our calculation differs most significantly from this earlier 
work in its treatment of the B ^ fo form factor, for which, for concreteness, we employ the CQM 
approach of Refs. [lOj HH 132], HSj. The earlier work simply asserts that the B ^ fo form factor 
ought be comparable to that of 5 — ttq. Our numerical studies support this, though we assign 
a large uncertainty to our assessment. The parameter scan technique we employ borrows heavily 
from the work of Beneke [H], in which the deviations to Sf{t) from sin(2/3) were studied for the 
B {7T^p'^,ri,ri',(f))Ks decay modes. That work eschewed the B foKg decay mode due to 
the uncertain qq structure of the /o [S3]; thus we have studied the possible numerical uncertainty 
incurred by this with great care. Our work also differs from Beneke's in that we explore the 
theoretical model space for input parameters which vary within 3a of their default values, as well 
as identify the parameter morphologies which give the largest excursions in AS^Ks- 

The assumed quark structure of the /o(980) does enter in the evaluation of the hadronic matrix 
elements and in the assessment of the /o decay constant and the S — > /o form factor in particular. 
We have investigated the range of ASf^Ks which results from varying either the B fo form factor 
and the /o scalar decay constant, after imposing the branching ratio constraint, giving [0.026, 0.028] 
for variations within la and [0.025,0.028] for variations within 3a. The value of ASfgXs simply 
not sensitive to the value of the F^-^" form factor and the /o decay constant. Although we can 
accommodate the experimental branching ratio in a qq model of the /o, we can not draw conclusions 
concerning the structure of the /o on the basis of our analysis. 

Let us summarize our results for the range of ASj^Ks ^foKsi which emerge from a scan of 
the complete theoretical parameter space, allowing for uncertainties in the infrared-cutoff depen- 
dent 0{Aqcd/Mb) corrections, which we characterize by the parameters Xa and Xh- In limited 
regions of the parameter space, the P'^ amplitude can become small, driving large excursions in 
ASff^Ks- Such excursions are removed, however, by demanding that the theoretical model confront 
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the empirical branching ratio up to some tolerance, that is, up to la for the Icr scans and up 
to 3(T for the 3a scans. Once we demand that our theoretical models satisfy the branching ratio 
constraint, the possible seed averaged range in ASf^Xs greatly reduced, giving [0.018, 0.033] and 
[—0.019, 0.064] for the scans in which the parameters are allowed to range within la and 3a of 
their central values, respectively. In comparison, we find that Cf^Ks ranges over [—0.013, —0.0024] 
and [—0.045,0.051], respectively. Generally, we find df^^^s to be sufficiently small that the values 
of ASf^Ks ^foKs are uncorrelated, note Fig. [ij It is worth emphasizing that we estimate the 
ranges using the extremal values of our parameter scans. This is a conservative approach. If we 
were to define the range of ASfgXs capture, e.g., 95% of the models about its most likely value 
in theory space, we would have much smaller ranges. Retaining our current approach, we could also 
sharpen our estimates once improved measurements of the branching ratio become available. 

Nevertheless, let us proceed to compare our range for ASj^Ks with the empirical result. The 
time-dependent CP asymmetry induced by the interference between BB mixing and direct decay, 
as well as the direct CP asymmetry, measured by Belle and BaBar [71 |Ml ES] are 0.85 ± 0.07 and 
0.08 ± 0.12, respectively. The value of sin(2/?) measured by Belle and BaBar [7] in i? ^ J /'^Ks 
decays and related charmonium modes is 0.668 ± 0.026, so that the empirical value of ASf^Xs 
is 0.18 ± 0.10 for the B —>■ f^Ks decay mode. Its error is substantially larger than our largest 
estimated value of ASf^Ksi so that further refinement of the experimental results is warranted. If 
we compare our results with those of Ref. [H] for the B — > {n^p^,ri,ri',(j))Ks decay modes, we see 
that the B fo^s mode compares favorably: it yields small values of AS*/, as do the B (pKs 
and B rj'Kg modes. Since both B — > foK^ and B (pKg decays lack color-suppressed tree 
contributions, the estimation of ASf should be particularly reliable, so that significant differences 
in the empirical values of ASf in these modes could truly signal new physics. 

We find that the largest single variation in our predictions of the decay amplitudes come from the 
uncertainty associated with the inverse moment of B meson distribution amplitude which enters 
the decay amplitude through the hard scattering terms. In our scans for which the parameters 
range over 3a of their central values, the hard scattering terms can dominate over all other terms in 
the decay amplitude and produce large values of the ^ amplitudes. The cancellation of a\ with 
Cg refiected in the default values of the parameter set is no longer possible for such large values 
of the hard scattering amplitude. This mechanism generalizes to any other mode where a\ and a\ 
interfere destructively. However, these large excursions in the A" ^ amplitudes do not generate large 
excursions in ASj^Ks^ they also act to enhance the charm penguin amplitude. Nevertheless, the 
extremal values of ASf^Ks^ negative values in particular, do come from the uncertainties 

associated with A^, Xa, and Xh- One expects this observation to be relevant to other b ^ s mode 
as well, though if there is a color-suppressed tree amplitude, there is more freedom for AS to be 
large. 

In future studies we would like to consider how the finite width of the /o(980) resonance can 
impact the value of ASf^Xs determined from B — >■ JqKs — > tt^tt'Ks and B foKs — > K^K~Ks 
decay processes. In particular, we wish to consider the role of the /o tt+tt^ and /o — * K^K^ 
scalar form factors, as well as that of possible OZI-violating effects, in the determination of ASj^Ks 
in a later publication [28] . 
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APPENDIX A: ANNIHILATION AMPLITUDES IN QCD FACTORIZATION 

The weak annihilation contributions 63 and bs^EW in the B foKs decay processes are ex- 
pressed as linear combinations of the annihilation amplitudes 3- Explicit expressions for these 
amplitudes in terms of the light-cone distribution amplitudes for scalar and pseudoscalar mesons are 
given in [TTl |29]. The general expressions for the leading-twist, light-cone distribution amplitudes 
for pseudoscalar and scalar mesons are written as 

{00 
n=l 

$5(x,/i) = /56x(l-x)|so + X] B,„(/i)C^/2(2a;_i)|. (Al) 

^ m=l ^ 

Here Bq is zero, as are all the even Gegenbauer moments, in the SU(3) limit. In order to present 
simple expressions for the A^{2 3 amplitudes, we truncate the expansion of the pseudoscalar meson 
after the second Gegenbauer polynomial and that of the scalar meson after the third Gegenbauer 
polynomial. 

We perform the integration over the light-cone distribution amplitudes in such a way that we 
always have I^y = lyx, where 



^xy = dx dyf{x,y) 



hy = I dy dxf{x,y). (A2) 







We retain finite e throughout the calculation and put e ^ only at the end to avoid dropping finite 
contributions. This amounts to a model of the power corrections which differs slightly from that 
employed previously [T71 SO], giving slightly different results for A\{PS) and Al{SP). Instead of 
the 7r^/3 found in Refs. [TTl ED], we get vr^/G. That is, the annihilation amplitudes can be written 
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as 



A[{PS) 



+53[(593 



fii[(180 - IStt^) + af (3726 - STStt^) + (27720 - 280871^)] 
- eOvr^) + af (37305 - 37807r2) + (714168 - 7236071^)] 



+(3Q;f - Saf + 3)[Si(18Xa - 36) + B^{mXA - 190)] - 2f^^r^Xl 



(A3) 



A[{SP) = TiasfKfl 5i[(-540 + 547r") + af(3726 - 3787r 



+53[(-5930 + 6007r^) + af (124350 - 126007r^) + (-1190280 + 1206007r^)] 



■(-13860 + 14047r2)] 



+(-3^3 - 3B,)[6{Xa - 1) + l8a^{XA - 2) + 12a^(3X^ - 8)] + 2f^r^Xl , (A4) 



K/ 



UlV-i 191 TT^ 

'/o ' X L-^"-^U^^A - ^^^A ~^>^ 60i?3(A^ - y^A + ^ + y )J 



^^(P5) = TiasfKfl <[18i?i(Xl- 4X^ + 4 



2 2 

+fJ[6(Xl - 2X^ + y ) - 18af (Xl - 4X^ + 4 + ^) 
+36af(Xl-^X^+^ + ^)]j 



(A5) 



TT 
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191 TT^ 



^^(5P) = Tra t"^ [-18^1 (XI - 4X^ + 4 + -) - 6053(^1 - + ^ + 



-ff [6(Xl - 2X^ + 



TT 



+ 18af (Xl-4X^ + 4+-) 



+36af(Xl-^X^+^ 



TT 



(A6) 



4(P5) = 7raJ;^/;„X^[rJ[6Pi(6X^-ll) + P3(120X^- 374)] 

-fJ[6(2X^ - 1) - eaf (6X^ - 11) + a^{12XA - 186)]] , 



(A7) 



Ai{SP) = 7rajK/XXA[rJ'[65i(6X^-ll) + S3(120X^- 374)] 
-fJ[6(2XA - 1) + 6Q;f (6Xa - 11) + af (72Xa - 186)] 



(A8) 



25 



[1] Y. Nir and H. R. Quinn, Ann. Rev. Nucl. Part. Sci. 42, 211 (1992); I. Dunietz, Annals Phys. 184, 

350 (1988) and references therein. 
[2] A.B. Carter and A.I. Sanda, Phys. Rev. Lett. 45 (1980) 952; Phys. Rev. D 23 (1981) 1567. 
[3] I.I. Bigi and A.I. Sanda, Nucl. Phys. B 193 (1981) 85; Nucl. Phys. B 281 (1987) 41; 
[4] I. Dunietz and J.L. Rosner, Phys. Rev. D 34 (1986) 1404. 

[5] For reviews, note Ref. [1^, as well as, e.g., J. Liu and L. Wolfenstein, Phys. Lett. B 197, 536 (1987); 
Y. Nir and D.J. Silverman, Nucl. Phys. B 345, 301 (1990); C. Dib, D. London, and Y. Nir, Int. J. 
Mod. Phys. A 6, 1253 (1991); T. Goto, N. Kitazawa, Y. Okada, and M. Tanaka, Phys. Rev. D 53, 
6662 (1996); N.G. Deshpande, B. Dutta, and S. Oh, Phys. Rev. Lett. 77, 4499 (1996); J.P. Silva and 
L. Wolfenstein, Phys. Rev. D 55, 5331 (1997); M.P. Worah, Phys. Rev. D 54, 2198 (1996); M. Gronau 
and D. London, Phys. Rev. D 55, 2845 (1997); Y. Grossman, G. Isidori, and M. P. Worah, Phys. Rev. 
D 58, 057504 (1998). 

[6] Y. Grossman and M. P. Worah, Phys. Lett. B 395, 241 (1997); R. Fleischer, Int. J. Mod. Phys. A 12, 
2459 (1997); M. Ciuchini, E. Franco, G. MartineUi, A. Masiero, and L. Silvestrini, Phys. Rev. Lett. 
79, 978 (1997); D. London and A. Soni, Phys. Lett. B 407, 61 (1997). 

[7] Heavy Flavor Averaging Group, |http://www.slac. stanford.edu/xorg/hfag . 

[8] G. Buchalla, G. Hiller, Y. Nir, and G. Raz, JHEP 0509,^0^4 (2005). 

[9] H. Boos, T. Manuel, and J. Renter, Phys. Rev. D 70, 036006 (2004). 
[10] H.-n. Li and S. Mishima, JHEP 0703, 009 (2007). 

[11] J. Charles et al. [CKMfitter Group], Eur. Phys. J. C 41, 1 (2005); note also 

|ht t p : / / www .slac.stanford.edu/xorg/ckmfitter|. 
[12] M. Bona et al. [UTfit Collaboration], JHEP 0507, 028 (2005); JHEP 0603, 080 (2006); JHEP 0610, 

081 (2006); note also 'http: / /www.utfig.org/ . 
[13] Y. Grossman, Z. Ligeti, Y. Nir, and H. Quinn, Phys. Rev. D 68, 015004 (2003). 
[14] M. Beneke, Phys. Lett. B 620, 143 (2005). 

[15] H. Y. Cheng, C. K. Chua, and A. Soni, Phys. Rev. D 72, 014006 (2005). 

[16] H. Y. Cheng and K. C. Yang, Phys. Rev. D 71, 054020 (2005). 

[17] H. Y. Cheng, C. K. Chua, and K. C. Yang, Phys. Rev. D 73, 014017 (2006). 

[18] R. L. Jaffe, Phys. Rev. D 15, 267 (1977). 

[19] J. D. Weinstein and N. Isgur, Phys. Rev. D 41, 2236 (1990). 

[20] D. Lohse, J. W. Durso, K. Holinde, and J. Speth, Phys. Lett. B 234, 235 (1990). 

[21] T. N. Truong, Phys. Rev. Lett. 61, 2526 (1988). 

[22] A. Dobado, M. J. Herrero, and T. N. Truong, Phys. Lett. B 235, 134 (1990). 
[23] J. Gasser and U.-G. Meii3ner, Phys. Lett. B 258, 219 (1991). 

[24] J. A. Oiler and E. Oset, Nucl. Phys. A 620, 438 (1997) [Erratum-ibid. A 652, 407 (1999)]. 

[25] U.-G. Meii3ner and J. A. Oiler, Nucl. Phys. A 679, 671 (2001). 

[26] T. A. Lahde and U.-G. Meii3ner, Phys. Rev. D 74, 034021 (2006). 

[27] S. Gardner and U.-G. Meii3ner, Phys. Rev. D 65, 094004 (2002). 

[28] R. Dutta, S. Gardner, T.-A. Lahde, and U.-G. Meifiner, in preparation. 

[29] M. Beneke, G. Buchalla, M. Neubert, and C. T. Sachrajda, Phys. Rev. Lett. 83, 1914 (1999); Nucl. 

Phys. B 606, 245 (2001). 
[30] M. Beneke and M. Neubert, Nucl. Phys. B 675, 333 (2003) and references therein. 



26 



[31] C. W. Bauer, S. Fleming, D. Pirjol, and I. W. Stewart, Phys. Rev. D 63, 114020 (2001). 

[32] C. W. Bauer, D. Pirjol, and I. W. Stewart, Phys. Rev. D 65, 054022 (2002). 

[33] C. W. Bauer, D. Pirjol, and I. W. Stewart, Phys. Rev. D 67, 071502 (2003). 

[34] G. Buchalla, A. J. Buras, and M. E. Lautenbacher, Rev. Mod. Phys. 68, 1125 (1996). 

[35] H.-n. Li, In the Proceedings of 5th Flavor Physics and CP Violation Conference (FPCP 2007), Bled, 

Slovenia, 12-16 May 2007, pp Oil [ arXiv:0707.1 294 [hep-ph]]. 
[36] C. W. Bauer, D. Pirjol, I. Z. Rothstein, and I. W. Stewart, Phys. Rev. D 70, 054015 (2004). 
[37] M. Beneke, G. Buchalla, M. Neubert, and C. T. Sachrajda, Phys. Rev. D 72, 098501 (2005). 
[38] A. Jain, I. Z. Rothstein, and I. W. Stewart, arXiv:0706.3399 [hep-ph]. 

[39] P. Colangelo, G. NarduUi, N. Paver, and Riazuddin, Z. Phys. C 45, 575 (1990); M. Ciuchini, E. Franco, 
G. Martinelh, and L. Silvestrini, Nucl. Phys. B 501, 271 (1997); M. Ciuchini, R. Contino, E. Franco, 
G. Martinelh, and L. Silvestrini, Nucl. Phys. B 512, 3 (1998) [Erratum-ibid.B531:656-660,1998]; 
M. Ciuchini, E. Franco, G. Martinelh, M. Pierini, and L. Silvestrini, Phys. Lett. B 515, 33 (2001); 
C. Isola, M. Ladisa, G. Nardulh, T. N. Pham, and P. Santorelh, Phys. Rev. D 64, 014029 (2001); 
S. J. Brodsky and S. Gardner, Phys. Rev. D 65, 054016 (2002). 

[40] R. Gatto, G. Nardulh, A. D. Polosa, and N. A. Tornqvist, Phys. Lett. B 494, 168 (2000). 

[41] A. Deandrea, R. Gatto, G. Nardulh, A. D. Polosa, and N. A. Tornqvist, Phys. Lett. B 502, 79 (2001). 

[42] A. Deandrea, N. Di Bartolomeo, R. Gatto, G. Nardulh, and A. D. Polosa, Phys. Rev. D 58, 034004 
(1998). 

[43] A. D. Polosa, Riv. Nuovo Cim. 23N11, 1 (2000). 

[44] D. Ebert, T. Feldmann, and H. Reinhardt, Phys. Lett. B 388, 154 (1996). 
[45] A. D. Polosa, private communication. 

[46] W. M. Yao et al. [Particle Data Group], J. Phys. G 33, 1 (2006). 

[47] Note Ref. HI]. We use the "Summer 2007" results. 

[48] A. Ah, B. D. Pecjak, and C. Greub, [MXiv:0709.4422 | [hep-ph] . 

[49] B. Aubert et al. [BABAR Cohaboration] , Phys. Rev. D 73, 031101 (2006). 

[50] A. Garmash et al. [Belle Collaboration], Phys. Rev. D 75, 012006 (2007). 

[51] D.E. Knuth, The Art of Computer Programming (TAOCP), Volume 2, Seminumerical Algorithms, 

Third Edition (Addison- Wesley, Reading, Massachusetts, 1997), p. 102. 
[52] M. Bona et al. [UTfit Collaboration], arXi v:0803.0659 [hep-ph]. 
[53] M. Beneke, private communication. 

[54] B. Aubert et al. [BABAR Cohaboration], arXiv:0708.2097 [hep-ex]. 
[55] K. Abe et al. [Belle Collaboration], Phys. Rev. D 76, 091103 (2007). 



27 



